clear; clc;
% clc;
%% Setup
InputParam_0    =   struct('XI',1/82.5/4,'Pir_SS',0.02/4,'KAPPA',0.00/4,...
                           'NU_C',1,'NU_L',1,...
                           'TrProb_ext',0.99,'Prob_ext',0.33,...
                           'TrProb_H',0.99,'Prob_H',0.37,...
                           'Fin_DomShare',0.79,'Fin_AdjCost',1,...
                           'Cons_FracT',0.5,'Cons_FracH',0.6,...
                           'B_Total_SS',0.91,'b_lb',-0.29,...
                           'WGap_ext',0,'WGap_H',0);

Param_Benchmark =   SubFun_SteadyState_OutputParam(InputParam_0);                       
%% Different Integration Integration 
% Financial
List_WGap_ext   =   (-0.25:0.05:0.25)';
Num_Fin         =   length(List_WGap_ext);
List_Param_Fin  =   cell(Num_Fin,1);
parfor ii=1:Num_Fin
    TempInputParam      =   InputParam_0;
    TempInputParam.WGap_ext=List_WGap_ext(ii);
    List_Param_Fin{ii}  =   SubFun_SteadyState_OutputParam(TempInputParam);
end

% Real 
List_WGap_H     =   (-0.25:0.05:0.25)';

Num_Real        =   length(List_WGap_H);
List_Param_Real =   cell(Num_Real,1);
parfor ii=1:Num_Real
    TempInputParam      =   InputParam_0;
    TempInputParam.WGap_H=  List_WGap_H(ii);
    List_Param_Real{ii} =   SubFun_SteadyState_OutputParam(TempInputParam);
    
end

List_PP         =   struct('Fin',{List_Param_Fin},'Real',{List_Param_Real});

save('Result_6_DifferentCorr.mat','Param_Benchmark','List_PP');


